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ABSTRACT 

Within the framework of the external shock model of gamma-ray bursts (GRBs) afterglows, we 
perform a morphological analysis of the early optical lightcurves to directly constrain model param¬ 
eters. We define four morphological types, i.e. the reverse shock dominated cases with/without the 
emergence of the forward shock peak (Type 1/ Type II), and the forward shock dominated cases with¬ 
out/with i/jn crossing the band (Type III/IV). We systematically investigate all the Swift GRBs that 
have optical detection earlier than 500 s and find 3/63 Type I bursts (4.8%), 12/63 Type II bursts 
(19.0%), 30/63 Type III bursts (47.6%), 8/63 Type IV bursts (12.7%) and 10/63 Type III/IV bursts 
(15.9%). We perform Monte Carlo simulations to constrain model parameters in order to reproduce 
the observations. We find that the favored value of the magnetic equipartition parameter in the for¬ 
ward shock {eq) ranges from 10“® to 10“^, and the reverse-to-forward ratio of eb {T^b) is about 100. 

The preferred electron equipartition parameter value is 0.01, which is smaller than the commonly 
assumed value, e.g., 0.1. This could mitigate the so- called “efficiency problem” for the internal shock 
model, if Eg during the prompt emission phase (in the internal shocks) is large (say, ~ 0.1). The 
preferred TZb value is in agreement with the results in previous works that indicates a moderately 
magnetized baryonic jet for GRBs. 

Subject headings: 


1. INTRODUCTION 

The first gamma-ray burst (GRB) afterglow emission 
was detected in 19 97, e.g. X-ray and optical afterglow 
from GRB 970228 (|Costa et aLllioOTtlvan Paradiis et al.l 
Il997fl . Over 18 years a variety of space- and ground- 
based facilities have detected hundreds of afterglows, 
with a wide coverage in both the spectral and tempo¬ 
ral domains (jKumar fc Zhanjl2015l for a recent review). 

The standard interpretation for the GRB afterglow 
emission was p roposed before the dis covery of the first 
afterglow data (iMeszaros fc Reeslll997ll . The general pic¬ 
ture is as follows ()Gao et al.ll2013l for a review): regard¬ 
less of the nature of progenitor and central engine, GRBs 
are believed to originate from a “fireball” moving at a 
relativistic speed. The fireball will inevitably be decel¬ 
erated through a pair of shocks (forward and reverse) 
propagating into the ambient medium and the fireball 
itself. Electrons are accelerated in both shocks and give 
rise to bright non-thermal emission through synchrotron 
or inverse Gompton radiation. Due to the deceleration of 
the fireball, a broad band afterglow emission with power- 
law rising and decaying behavior is expected for the GRB 
afterglow. 

In the pve-Swift era, the simple external shock 
model provided successfu l interpretations for a large 
array of afterglow data (IWiier s et al.l I l997t IWaxmad 
199’^ IWiiers fc Galamal II999I: iHuang et al.lll999l 120001 : 

Panaitescu fc Kumaill200lL 120021 : lYost et al.ll2003ll . al- 
though moderate revisions were sometimes required, for 
instance, invoking w ind-type density medium instead 
of constant density (|Dai fc Lul 1199811 IMeszaros et al.l 


Il998t iGhevalier fc LilfTgM . 1200011 . refi ning the joint for¬ 
ward s hock and reverse shock s ignal (IM eszaros fc ReesI 
I997t IM eszaros fc R.eesI II999I: ISa ri fc Pirad Il999allbl: 


Kobavashi fc Zhai^ 120031: IZhang et al.l 120031: IWu et all 
20031: IZou et al.l 1200511 . considerin g continuous en- 


2UUdt IZou et al.l l2UUi3l) . considerin g: continuous en¬ 
ergy injection into the blastwave IPai fc Liil Il998al : 
iRees fc MeszarosllI998l:[^ang fc Meszarosll200Ill. takin; 


_. - _... __ _^ _._ _ _g' 

into account the jet break effect (iRhoadd I999HSari et al-l 

II999I: IZhang fc MeszarosI 120021 : iRossi et al.l 1200211 . etc. 
Entering the Swift era, some new unexpected signatures 
in GRB aftergj^ l ows w e re revealed (iTagliaferri et al.ll2005l: 
Burrows et al.! 1^00^ IZhang et al.l 120061: Nousek et al l 


2006HO’Brien et al.ll2006l:' Evans et al .112009 ), which how¬ 

ever are still acommodated within the standard frame¬ 
work, provided some additional physical processes are 
invoked self-consistently, such as a late centr al engine 
activity (jZhang et al.ll2006l: iNoiisek et al.ll2006ll . 

The external shock model is elegant in its simplicity, 
since it invokes a limited number of model parameters 
(e.g. the total energy of the system, the ambient density 
and its profile), and has well defined predicted spectral 
and temporal properties. Given this model, the accu¬ 
mulation of afterglow data has led to great advances in 
revealing physical properties in GRB ejecta as well as 
the circum-burst medium. In practice, there are two 
approaches for applying the external shock model to 
the observational data: one can start with the data, fit 
the lightcurve and spectrum with some empirical broken 
power-law functions to get both a temporal index a and a 
spectral index /3, and then constrain the related afterglow 
parameters by applying the so called “closure relation” 
(|Zhang fc Mdszaro4l2004 iGao et all 120131 : IWang et al.l 
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HqH); or alternatively, one can start with the external 
shock model, draw predicted lightcurves and spectrum 
with varying parameters, and then constrain the rele¬ 
vant parameters by fitting the observational data with 
the theoretical prediction. 

Both approaches encounter their own difficulties. The 
first approach is usually non-optional, since some compli¬ 
cated effects such as the equal arrival time effect and the 
gradual evolution of cooling frequency res ult in a smooth¬ 
ing o f the spectral and tem poral breaks (iGranot &: Saril 
I 2 OO 2 I : lUhm fc Zhand l2014f) , leading to imperfection of 
the “closure relation”. For the latter approach, due to 
the simple behavior of the afterglow data and the sim¬ 
ple power-law property of the synchrotron external shock 
model, the model parameters obtained by fitting indi¬ 
vidual bursts usually suffer severe degeneracy, unless the 
observed SED could fully cover all synchrotron charac¬ 
teristic frequencies (Kumar & Zhang 2015 for a review), 
e.g. Va (self-absorption frequency), Vm (the characteristic 
synchrotron frequency of the electrons at the minimum 
injection energy), and Vc (the cooling frequency). For the 
cases that all the observations are in the same spectral 
regime or only covers one break frequency, which usually 
happened when only optical and X-ray data are avail¬ 
able, individual fitting cannot make tight constraint on 
the parameters such as Ce and cb, even though the model 
calculated lightcurve could nicely fit the data (e.g. see 
the recent results presented in Japelj et al. 2014). 

Regardless of these difficulties, both approaches work 
best with the observational data in the early stages, 
which contain much richer information. On the other 
hand, although multi-wavelength observations are rou¬ 
tinely carried out for many GRBs, optical data are still 
generally the most valuable for model constraint. The 
reason is as follows: first, data in the radio band are still 
limited, although several radio lightcurves have been in¬ 
terpreted in detail and s ome interesting studies h ave been 
performed statistically (IChandra fc Frail I2012L and ref¬ 
erence therein). Second, the X-ray lightcurves are some¬ 
times dominated by the X-ray flares and X-ray plateaus, 
which can not be fully interpreted with the simple ex¬ 
ternal shock model, additional physical processes, such 
as a radiation component that is r elated to the late 
central engiiie activity, are needed (IZhang et al.l 120061: 
Nousek et al.l l2006t iGhisellini et al.l [2007t Kumar et all 
2008bl laf). Last but not least, in the early stage, the re¬ 
verse shock spectrum is expected to peak in the optical 
band. Investigation of reverse shock emission is very im¬ 
portant for studying the detailed features of GRB ejecta, 
such as the composition of the jet, since its radiation 
comes directly from the shocked ejecta materials. 

Assuming that afterglow parameters for different 
GRBs come from the same distributions, statistical prop¬ 
erties of a sample of GRBs can be used to constrain the 
global features of model parameters . For instance, in 
the pre-Swift era, IZhang et al.l (l2003ll proposed to cate¬ 
gorize the early optical afterglows into different types of 
combinations of reverse and forward shock emission and 
they suggested that the afterglow parameter space could 
be explored based on a morphological analysis. After a 
decade of successful operation of Swift , a fairly good 
sample of early afterglow lightcurves is in hand. It is 
now of great interest to develop and implement the mor¬ 
phological analysis on the current observations to make 


reasonable constraints on the model parameters. Specif¬ 
ically, the morphological analysis method can be divided 
into two separate parts: Monte Garlo simulations and ob¬ 
servational sample analysis. In the simulation part, we 
assume some intrinsic distributions for each afterglow pa¬ 
rameter, simulate a sample of afterglow lightcurves, and 
then distribute them into their relevant lightcurve types. 
In the sample analysis part, we try to find a well defined 
sample of early optical lightcurves, and calculate the rel¬ 
ative number ratios among different lightcurve types. By 
comparing the results between these two parts, one can 
make constraints on relevant parameters. 

The structure of the paper is as follows. We illustrate 
the morphological analysis method and the sample selec¬ 
tion process in section 2, including the definition of differ¬ 
ent lightcurve categories, and the theoretical scheme for 
determining categories for given values of the afterglow 
parameters. In section 3, we apply the morphological 
analysis to the GRB sample with Monte Garlo simula¬ 
tions, and explore the parameter regimes by comparing 
the simulation results and observations. We discuss our 
results in section 4, and briefly summarize our conclu¬ 
sions in section 5. Throughout the paper, the convention 
Q = lO^Qn is adopted for cgs units. 

2. EARLY OPTICAL AFTERGLOW MORPHOLOGY 
2.1. Lightcurve classification 

The morphology of early optical afterglows essentially 
reflects the relative relation between the forward shock 
and the reverse shock emission. Since the strength of the 
forward and reverse shocks are determined by the same 
set of GRB parameters, namely the initial Lorentz factor, 
the kinetic energy of the fireball, the circum-burst density 
and the microphysics parameters, in principle a study of 
the morphology can yield direct model constraints. 

In previous works, the early optical afterglows for con¬ 
stant density me dium model were usually classified into 
three categories (IZhang et al.ll2003L iJin fc FaDll2007D : 

• Type I: re-brightening. At the very early stage, the 
lightcurve is dominated by the reverse shock emis¬ 
sion, but later a re-brightening signature emerges 
due to the forward shock emission contribution. 
Both reverse shock peak and forward shock peak 
are evident in this type of lightcurve. 

• Type II: flattening. The forward shock peak is be¬ 
neath the reverse shock component. The forward 
shock emission only shows its decaying part at the 
late stage, since the reverse shock component fades 
more rapidly. 

• Type III: no reverse shock component. In this 
case, the reverse shock component is either too 
weak compared with the forward shock emis¬ 
sion, or is completely suppressed for some rea- 
son, such as magnetic fields dominating the ej ect a 
(|Zhang fc Kobavashill2005l : iMimica et al.ll201^ . 

Note that for forward shock dominated cases (Type 
III), there are still two distinct shapes of lightcurve, de¬ 
pending on whether i/^(tx) is larger tha n r'nnt, or not, 
where tx is the reverse shock crossing time (|Sari fc Piranl 
dH). (t X) > i^opt, the rising slope of the lightcurve 
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(a) (b) 



Fig. 1.— Example light curves for all four types with typical afterglow parameters. Type h E = 10®^ erg, Fq = 100, u = 10 cm 
= 0.1, = 10-'^, 7^B = 100 and p = 2.3; Type II: E = lO®^ erg, Fq = 100, n = 0.01 cm"®, = 0.01, = 10"®, TZb = 100 

and p = 2.1; Type III: E = 10®^ erg, Fq = 100, n = 10 cm“®, = 0.01, eg = 10“^, TZb = 100 and p = 2.3; Type IV: E = 10®^ erg, 

Fq = 100, n = 10 cm“®, e^’^ = 0.1, eg = 10“®, TZb = 1 a-nd p = 2.3. 


would have a very clear steep or to shallow 

transition, otherwise the rising slope is always 
steep. Since an insight on the t'm(tx) value could lead to 
strong constraints on relevant afterglow parameters, in 
this work we further categorize the forward shock domi¬ 
nated lightcurves into two categories: 

• Type III: forward shock dominated lightcurves 
without I'm crossing. The observed optical peak 
is the deceleration peak. 

• Type IV: forward shock dominated lightcurves with 
I'm crossing. The observed optical peak is the Vm 
crossing peak. 

Figure [T] shows the example light curves for all four types 
with typical afterglow parameters. 

2.2. Theoretical scheme for determining categories 

Consider a uniform relativistic shell (fireball ejecta) 
with an isotropic equivalent energy E, initial Lorentz fac¬ 
tor Tq, and observed width Aq, expanding into a homo¬ 
geneous interstellar medium (ISM) with particle number 


density n at a redshift z. During the initial interaction, a 
pair of shocks develop: a forward shock propagating into 
the medium and a reverse shock propagating into the 
shell. After the reverse shock crosses the shell (at tx), the 
forward shock follows the Blandfor d-McKee self-similar 
solution (|Blandford fc McKedllQT^ . Synchrotron emis¬ 
sion is expected behind both shocks, since electrons are 
accelerated at the shock fronts via the Ist-order Fermi ac¬ 
celeration mechanism, and magnetic fields are believed 
to be generated behind the s hocks due to plasma in- 
stabilities (for forward shock) (iMedvedev fc LoebllI999ll 
or shock compression amplification of the magnetic field 
carried by the central engine (for reverse shock). The in¬ 
stantaneous synchrotron spectrum at a given epoch can 
be described with three characteristic frequencies t'a, , 
and Vc, and the peak synchrotron flux density Fb^max 
(jSari et al.lll998h . The evolution of these four parame¬ 
ters can be calculated, with the help of notations to pa¬ 
rameterize the microscopic processes, i.e. the fractions 
of shock energy that go to electrons and magnetic fields 
(ce and cb), and the electron spectral index p. It is then 
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straightforward to calculate the lightcurve for a given 
observed frequency (e.g. optical frequency t'opt), 


,E,ro,Ao), (1) 


where the superscript r and / represent reverse and for¬ 
ward shock respectively. 

In principle, the morphology for a specific GRB can 
be determined once the entire optical lightcurve is calcu¬ 
lated. However, this process is time consuming and not 
conducive for realizing Monte Carlo simulation to ex¬ 
plore a large parameter space. However, exploiting the 
power-law behavior of the afterglow emission, here we 
propose an efficient scheme to categorize the lightcurve 
type, by comparing the reverse shock and forward shock 
flux strength at some special time point, rather than 
comparing them for the entire duration. The detailed 
scheme is illustrated as follows: 

The dynamical evolution during the rev erse shock 
cross ing phase can be classified into two cases (|Kobavashil 
l200(tl . depending on whether the reverse shock becomes 
relativistic in the frame of the unshocked shell mate¬ 
rial (thick shell case) or not (thin shell case). Since the 
emissions from both reverse shock and forward shock be¬ 
have differently in each case, our first step is to deter¬ 
mine the thin/thick properties for given set of parame¬ 
ters. One practical way is to compare the duration of 
the burst T = Aq/c and the deceleration time of the 
ejecta tdec = (3i?/327rnmprQC®)^/^, i.e., T > tdec is for 
thick shell case an d T < is for t hin shell case, and 
tx = max(tdec, T’) (IZhang et al.ll2003h . 

For the thin shell case, before shock cros sing time, 
the evolution of , and T'J'.max reads^ (|Cao et al.l 

l2nT3h 


< = 1.9 X 10^^ 
= 4.1 X 10^® 


Hz z 


^52 ^ 0 ( 2 ^ 


-7 p-2t^18 

G(2.3) 

Hz 

nz^i0,2"'0,0 ^B,-2''2 


5/2 2 

0.0 


. /2 . 
B,-2 


6 

2i 


F'' 

t',max 


9.1 X 10® /rJy z ^/^£^ 52 ^r^, 2 »^o.oeB .-2 


D 


- 2 , 3/2 
28 ^2 


yf _ Q 1 V inl6 u, --l_^(p)_p4 „l/2 2 1/2 

X nz Z 0,2^0,0 ^e,-l^B,-2y 


izi = 4.1 X 10®® Hz 


'1 = 

^/max = 1-1 X 10^ /rJy z 


8 „3/2 1/2 ^-2,3 

0,2^^0,0 -2-^28 ^2’ 


( 2 ) 

where G{p) = = (1 -|- z) is the redshift 

correction factor. For simplicity, we omit the superscript 
of and e^-^. With the expression of the values of 
and at the shock crossing time could 

be easily obtained. Consequently, with the te mporal evo¬ 
lution power-law indices of these parameters (|Cao et al.l 
1201,311 . one can calculate their values for the post shock 
crossing phase. At txj we have 


) _ 2 c; V in5 "8^-2/3p-2/3 -4/3 -2 -2 

,,r (4 \ ~ X iU Z £,52 4 0,2 ^0.0 ^e,-l^B.-2> 

f ,, . — zo.o Z ^52 4 0.2 "-0.0 2- 

Vin(tx) 


^ Since we focus on afterglow emission in optical band, the effect 
of Ua is not considered here. 


(3) 

where G{p) factor is normalized to p = 2.3. We can see 
that for the time we are interested in (e.g., mainly around 
or after tx), both reverse shock and forward shock emis¬ 
sion would be in the “slow cooling ” regime (i^g > for 
reasonable parameter regimes ^ (|Sari et al.lll998D . We 
take slow cooling for both reverse and forward shock 
emission in the following, so that the shape of the light 
curve essentially depends on the relation between 
and r'opt (similar arguments also apply to the thick shell 
case). The evolution of for the thin shell case is 
shown in Figure [2^, which reads 

(Xt° {t< tx), 4-1 cx {t > tx), 

oc t® (t < tx), v'm fx {t > tx)- (4) 

When 4'^(t X) is larger than r-opt, we call it FS I case (oth¬ 
erwise FS H case), and would cross the optical band 
once (at tf). Similarly when v^{tx) is larger than r'opt, 
we call it RS I case (otherwise RS H case), and r-Jjj would 
cross the optical band twice (at C,i arid tr, 2 )- There are 
altogether four combinations for different shapes of re¬ 
verse and forward shock light curves. For each combi¬ 
nation, we first check if the peak of the reverse shock 
emission is suppressed by the forward shock. If so, the 
lightcurve belongs to Type HI (FS H) or IV (FS I). Oth¬ 
erwise, we will further check if the peak of the forward 
shock emission is suppressed by the reverse shock. If so, 
the lightcurve belongs to Type H, otherwise it is Type I. 
The evolution of F^'f for all four cases are presented in 
Table 1 and shown in Figure [5)3. The scheme to catego¬ 
rize the lightcurve type is presented in Table 2. 

For the thick shell case, before shock crossing time tx, 
the evolution of , and F^’^^ reads 


= 7.6 X 10®® Hz z 


11 tr„ 11-1 t-2 „1/2 2 1/2 


G(2.3)^0’2"'0’° ee.-ieB.-2> 


i/g — 1.2 X 10®®" HZA 52 ^ ^o{l3.'^Ofl^B,-2^2 ^ 

= 1-3 X 10® pJy 

ui = 1.0 X 10®® 

4 = 1.2 X 10®" 


— 1-2 X 10® /iJy zA 52 ^ 0 , 13 ’^ofo ^b(- 2-^28 • 


(5) 


The evolution of for the thick shell case is shown in 
Figure [2J:, which reads 


4-^ oc t ^ {t <tx), vl, fxt {t > tx), 

(Xt° {t < tx), 4-;, OC i-'^3/48 ^ 


Similar to the thin shell regime, we define four cases for 
different evolution, and present the results for all 
cases in Table 1 and Figure [2Ji. The scheme to catego¬ 
rize the lightcurve type is presented in Table 2. Note that 


^ Note that for some extreme parameters, the fast cooling regime 
(iZc < I'm) might be relevant at the shock crossing time. However, 
in those cases, Uc{tx) could not be much smaller than i'm(tx)^ so 
that the real lightcurve shape would not deviate too much from 
the ones presented here (under the slow cooling assumption). 
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Thick shell case 



(c) 


Thick shell case 



(d) 


Fig. 2.— Illustration of the Vrn evolution (left panels) and optical lightcurves (right panels) for both forward shock (blue lines) and 
reverse shock (black lines) emission. Top panels are for thin shell regime and bottom panels are for thick shell regime. For Um evolution 
panels, red solid line represents the observer frequency. For all panels, solid lines are for cases with t'm^(fx) > I'opt, dot-dash lines are for 
cases with (tx ) < i^opt- Red circles on lightcurve indicate the points for comparison in order to categorize the lightcurve types. Dotted 
lines at the end of reverse shock indicate the high latitude emission after izj crosses the observer frequency. 


TABLE 1 

The evolution of for all cases 


Thin Shell 

Thick Shell 

FS I case (tf > tx) 

FS I case (tf > tx) 


t < tx 

tx < t < tf 

t > tf 

t < tx 

tx < t <tf 

t > tf 

fI oc 

t3 

tl/2 

t-3(p-l)/4 

H/3 

tl/2 

*-3(p-l)/4 

FS II case (tf < tx) 

FS II case (tf < tx) 


t < tx 

t > tx 

t < tf 2 

tf.2 <t <tx 

t > tx 

fI oc 

t3 

f-3(p-l)/4 

*4/3 

*(3-p)/2 

*-3(p-l)/4 

RS I case (tr > tx) 

RS I case (tr > tx) 


t < tr.l 

tr-,1 < t < tx 

tx t <C tj-^^ 

t > tr,2 

t < tx 

t X t <C. tr 

t > tr 

F[,^ 

j(6p-3)/2 

t-l/2 

f-16/35 

f-(27p+7)/35 

tl/2 

*-17/36 

*-(73p+21)/96 

RS II case > 

tr tx) 

RS II case (tr < tx) 


t < tx 

t>tx 

t < tx 

t > tx 

K oc 

t(6p-3)/2 

f-(27p-l-7)/35 

tl/2 

*-(73p+21)/96 


for thick shell case, Type III lightcurve could mimic like Type IV when t/ 2 ^ tx- However, the parameter space 
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TABLE 2 

The scheme to categorize the lightcurve type for all combinations for different shapes of reverse and forward shock 

LIGHT CURVES 


Thin Shell 

Thick Shell 

FS I + RS I 

FS I + RS I 


FHitr.l) > Fl{tr;l) 

f;Px) < Fl(t^) 

F7Px) > Fjpx) 

Type IV 

Flitf) > FKtf) 

Fi(ti) < FKtj) 

Type IV 



Type I 

Type II 

Type 1 

Type II 

FS I -F RS II 

FS I -F RS II 


E;(tx) > Fl(tx) 

F7(ix) < Fl(tx) 

F7(ix) > Flpx) 

Type IV 

Fl(tj) > F:(tf) 

Fi(tf) < F:(tj) 

Type IV 

Fi(tf,i) > F;(tj,i) 


Type I 

Type II 

Type I 

Type II 

FS II-F RS I 

FS II -F RS I 


F:<,tr,l) > Fi(tr,l) 

F7Px) < Fj(tx) 

f:px) > Fjpx) 

Type III 

Fi(tv,) > 

X 

A 

X 

Type III 

F,l(tx) > F;(tx) 

F,l(tx) < F;(tx) 

Type I 

Type II 

Type I 

Type II 

FS II + RS II 

FS II + RS II 

E;(tx) < Fj(tx) 

f;Px) > Flpx) 

F;(tx) < Fipx) 

F;(tx) > Fj(tx) 

Type III 

Fjpx) > F;(tx) 

A 

X 

A 

X 

Type III 

FxVtx) > F;(tx) 

Fj(tx) < F;(tx) 

Type I 

Type II 

Type I 

Type II 


for this situation is very limited and practically this con¬ 
fusion could be easily clarified with spectral information, 
thus we still count this case as Type III in the simulation 
results. 


2.3. Sample Selection 

We systematically investigate all the Swift GRBs that 
have optical detections at times earlier than 500 s after 
the prompt emission trigger, from the launch of Swift 
to March 2014. A sample of 114 lightcurves is compiled 
either from published papers or from GCN Circulars if no 
published paper is available (|Li et al.ll2012l : iLiang et"^ 
iMllKann ei~aIIIM0l[20nll . 

We first find the bursts without a detected initial ris¬ 
ing. We fit their initial decaying phase with a single 
power-law function, and keep the bursts which have a 
decaying slope larger than 1.5 as candidates for Type I 
or Type II. Other bursts with relatively slower slopes are 
excluded from the following analysis since in principle 
they could belong to any one of the four types. 

Within the remaining sample, we find the bursts with 
rebrightening or flattening (steep decay to shallow decay) 
features. For these bursts, we fit their initial rising and 
decaying part with a smooth broken power-law function 
and take the bursts with decaying slope larger than 1.5 
as candidates for Type I or Type II. All other bursts are 
taken as the candidate for Type III or Type IV. 

For Type I/II candidates, we fit their lightcurves with 
two separate broken power-law components. If the peak 
flux for the weaker component is completely suppressed 
by the stronger component, the lightcurve is classified as 
Type II, otherwise it is classified as Type I. For Type 
III/IV candidates, we fit their lightcurves with one bro¬ 
ken power-law component, and take a rising slope smaller 
than 0.6 as the division between Type III and Type IV. 
For some bursts, there are early observations that can 
be used to exclude Type I and Type II, but it is hard to 
determine their rising slope to justify a Type III or Type 
IV classification in some cases, for instance, if the data 
points are too close to the peak or if the rising phase 
is superposed on an optical flare. We count these as an 
overlapping type in the following analysis. 


Eventually, we find 3 Type I bursts (4.8%), 12 Type II 
bursts (19.0%), 30 Type III bursts (47.6%), 8 Type IV 
bursts (12.7%) and 10 Type III/IV bursts (15.9%). The 
lightcurves for each type are shown in Figure [3] and their 
properties are collected in Table 3, including the GRB 
name, the onset rising slope, decaying slope, peak time, 
peak flux and lightcurve type. 


3. MONTE CARLO SIMULATIONS 

In principle, if the intrinsic distribution function for 
each parameter listed in Equation [T] is available, we can 
simulate a sample of afterglow lightcurves, and distribute 
them into their relevant categories with the aforemen¬ 
tioned theoretical scheme. The properties of the parame¬ 
ter distributions could in turn be constrained by compar¬ 
ing the simulation results with the observational results 
collected in section [2^ 

In the literature, several statistical works have been 
done through fitting indivi d ual bursts, either for late 
(iPanaitescu fc Kumaill200lL l2002t lYost et al.l l2003ll or 


early ( Liang et al.l 120131 : iJaneli et al.ll2014D broad-band 
observations. Although the intrinsic distribution func¬ 
tions are still poorly understood, some useful informa¬ 
tion, such as typical values or distribution ranges for 
most of the afterglow parameters, have been proposed 
(jKumar fc Zhansl2015fl . With this information, we could 
assume some proper distributions for the afterglow pa¬ 
rameters, e.g. Gaussian distribution around the typical 
value or uniform distribution within a reasonable range. 
Even if the adopted distribution functions may be de¬ 
viated from the intrinsic ones, it is still possible to jus¬ 
tify how each parameter affects the result of morpholog¬ 
ical analysis. Nevertheless, for critical parameters that 
severely affect the results, their preferred values could be 
explored by comparing with the current observations. 


3.1. Simulation Setup 

The adopted distributions for generating afterglow pa¬ 
rameters listed in Equation [T] are as follows: 


• The redshift z is generated based on the assump¬ 
tion that the GRB rate roughly traces the star for- 
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TABLE 3 

Early optical afterglow properties for GRBs in the selected sample. 


GRB a“ 






/(.(xlO-^)" 

Typ' 

021004 1.50 ± 

0.16 

1.06 

± 

0.11 

100 

2.58 

I 

050525A 1.50 ± 

0.20 

1.15 

± 

0.10 

66 

23.23 

I 

090424 1.51 ± 

0.25 

0.85 

± 

0.15 

176 

0.35 

I 

990123 2.52 ± 

0.42 

1.52 

± 

0.15 

42 ± 9 

698.64 ± 140.00 

II 

021121 1.97 ± 

0.15 

1.05 

± 

0.11 

130 

0.20 

II 

050904 3.00 ± 

0.27 

1.20 

± 

0.12 

405 ± 40 

57.90 ± 15.40 

II 

060111B 2.41 ± 

0.21 

1.19 

± 

0.11 

29 

6.16 

II 

060117 2.42 ± 

0.11 

1.00 

± 

0.09 

109 

168.91 

II 

060908 1.48 ± 

0.25 

1.05 

± 

0.09 

50 

3.20 

II 

061126 2.00 ± 

0.32 

0.86 

± 

0.11 

23 

31.44 

II 

080319B 2.74 ± 

0.42 

1.23 

± 

0.15 

33 ± 12 

18246.50 ± 7200.00 

II 

081007 1.90 ± 

0.21 

0.68 

± 

0.08 

3 ± 1 

0.68 ± 0.16 

II 

090102 1.84 ± 

0.18 

1.10 

± 

0.11 

50 ± 11 

8.07 ± 3.00 

II 

091024 2.20 ± 

0.21 

1.05 

± 

0.11 

430 ± 52 

3.50 ± 1.10 

II 

130427A 1.67 ± 

0.11 

1.01 

± 

0.07 

13 ± 3 

2359.09 ± 500.00 

II 

030418 -0.81 ± 0.12 

0.55 

± 

0.04 

1190 ± 109 

0.35 ± 0.01 

III 

050820A -2.00 ± 0.21 

1.10 

± 

0.18 

500 ± 30 

2.00 ± 0.30 

III 

060110 -1.20 

0.80 

50 

15.00 

III 

060210 -1.19 ± 

0.18 

1.33 

± 

0.02 

718 ± 23 

0.10 ± 0.01 

III 

060418 -1.80 ± 

0.10 

1.20 

± 

0.03 

158 ± 1 

9.22 ± 0.09 

III 

060607A -2.39 ± 

0.18 

1.31 

± 

0.11 

180 ± 15 

4.50 ± 0.12 

III 

060926 -4.02 ± 

1.75 

0.75 

± 

0.12 

78 ± 10 

0.17 ± 0.03 

III 

061007 -2.99 ± 

0.03 

1.67 

± 

0.02 

90 ± 2 

179.09 ± 0.86 

III 

070419A -1.00 ± 

0.12 

1.28 

± 

0.03 

753 ± 32 

0.06 ± 0.01 

III 

070420 -1.43 ± 

0.54 

0.90 

± 

0.08 

196 ± 22 

1.80 ± 0.14 

III 

071010A -1.06 ± 

0.07 

0.72 

± 

0.01 

384 ± 30 

0.46 ± 0.02 

III 

071010B -0.70 ± 

0.37 

0.52 

± 

0.03 

147 ± 36 

0.32 ± 0.02 

III 

071025 -1.20 ± 

0.18 

1.11 

± 

0.12 

563 ± 80 

0.07 ± 0.01 

III 

071031 -0.74 ± 

0.02 

0.76 

± 

0.03 

1055 ± 11 

0.10 ± 0.01 

III 

080319A -1.80 ± 

0.08 

0.65 

± 

0.07 

238 ± 17 

0.02 ± 0.01 

III 

080603A -3.85 ± 

0.13 

1.17 

± 

0.02 

482 ± 14 

1.04 ± 0.05 

III 

080710 -1.20 ± 

0.12 

0.65 

± 

0.07 

1695 ± 42 

0.47 ± 1.10 

III 

080810 -1.15 ± 

0.05 

1.14 

± 

0.03 

142 ± 2 

14.28 ± 0.21 

III 

080928 -0.78 ± 

0.05 

2.18 

± 

0.19 

3223 ± 130 

0.40 ± 0.01 

III 

081008 -2.20 ± 

0.09 

1.09 

± 

0.04 

175 ± 1 

8.18 ± 0.06 

III 

081126 -1.14 ± 

0.02 

0.39 

± 

0.01 

159 ± 2 

1.55 ± 0.01 

III 

081203A -0.96 ± 

0.03 

1.37 

± 

0.02 

376 ± 1 

14.35 ± 0.03 

III 

090313 -1.36 ± 

0.19 

0.92 

± 

0.02 

1002 ± 69 

0.70 ± 0.04 

III 

090812 -1.35 ± 

0.32 

1.37 

± 

0.29 

71 ± 8 

1.79 ± 0.11 

III 

091029 -3.10 ± 

0.23 

0.48 

± 

0.07 

312 ± 52 

0.14 ± 1.00 

III 

100219A -1.50 ± 

0.54 

0.95 

± 

0.01 

619 ± 76 

0.15 ± 0.02 

III 

100901A -1.87 ± 

0.13 

1.00 

± 

0.07 

1200 ± 95 

0.18 ± 0.01 

III 

100906A -1.76 ± 

0.04 

1.10 

± 

0.03 

137 ± 1 

17.65 ± 0.17 

III 

110213A -1.92 ± 

0.02 

0.73 

± 

0.04 

230 ± 1 

4.64 ± 0.01 

III 

121217A -1.80 ± 

0.12 

0.80 

± 

0.01 

1806 ± 29 

0.03 ± 0.01 

III 

060218 -0.35 ± 

0.01 

0.78 

± 

0.03 

55032 ± 1305 

0.03 ± 0.01 

IV 

060605 -0.47 ± 

0.05 

1.24 

± 

0.01 

701 ± 38 

1.60 ± 0.05 

IV 

070318 -0.59 ± 

0.11 

1.07 

± 

0.04 

456 ± 25 

1.64 ± 0.04 

IV 

070411 -0.58 ± 

0.14 

0.96 

± 

0.01 

655 ± 25 

0.18 ± 0.01 

IV 

080330 -0.25 ± 

0.05 

0.97 

± 

0.03 

902 ± 45 

0.21 ± 0.01 

IV 

090510 -0.53 ± 

0.17 

1.13 

± 

0.21 

1273 ± 501 

0.07 ± 0.01 

IV 

120815A -0.25 ± 

0.03 

0.64 

± 

0.01 

521 ± 21 

0.09 ± 0.01 

IV 

120119A -0.28 ± 

0.08 

2.38 

± 

0.90 

1021 ± 63 

0.27 ± 0.01 

IV 

060729 -1.20 ± 

0.21 

0.90 

± 

0.15 

800 ± 82 

0.60 ± 0.02 

h 

060904B -0.95 ± 

0.11 

1.07 

± 

0.07 

493 ± 26 

0.40 ± 0.01 

h 

060906 -0.20 ± 

0.45 

1.03 

± 

0.35 

1263 ± 639 

0.05 ± 0.01 

h 

070611 -2.58 ± 

0.56 

0.90 

± 

0.21 

2114 ± 342 

0.09 ± 0.01 

h 

071112C -0.60 ± 

0.37 

0.91 

± 

0.02 

165 ± 13 

0.32 ± 0.02 

h 

080319C -0.38 ± 

0.05 

2.15 

± 

0.10 

654 ± 39 

0.21 ± 0.01 

h 

081109A -0.19 ± 

0.18 

0.94 

± 

0.03 

559 ± 128 

0.25 ± 0.03 

h 

090726 -1.27 ± 

0.12 

0.70 

± 

0.18 

290 ± 45 

0.12 ± 0.03 

h 

110205A -3.54 ± 

0.42 

1.51 

± 

0.15 

958 ± 56 

3.24 ± 0.33 

h 

120711A -0.50 ± 

0.10 

1.03 

± 

0.12 

332 ± 2 

7.15 ± 1.05 

h 


Data Reference 
^iraba^^^L ^003) 

(^20^) 

^astro^irad^^^^^ (^^9^ 
(^200^ 

{^00^ 

^tratta^^^L ^200^ 
^elin^^^^^L ^006) 
(^20^) 
(^200^ 
^200^ 

^0^) 

^^estran^^^^L ( 20 ^^ 
^00^ 
^200^ 
lA (2006) 

^urrar^^^L (^00^ 
^olinar^^^L ^00^ 
^olinar^^^L ^00^ 
(^200^ 

^^indel^^^L (^00^ 
^elandr^^^j^ ^200^ 
(^2008) 


^008) 

^200^ 

^erie^^^ab 
Kriihler et al. ('2009a) 



Cenko (2008) 

^Addorz^^^l^ (2011) 

Kruhle^e^^ 

Ross^^aL 

^Qi^^e^aL 

Klot^e^aL 

Kuh^^^. 

Melandr^^^ 

^00^ 

^arshal^^^ru^^ ^00^ 
(^20^) 

^orbov^ko^^^^l^ (2G1^ 
^orbovsko^^^^l^ (2G1^ 
^ucchiar^^^^j^ ^0^^) 
^lliot^^^^L ^OM) 
^200^ 
{^00^ 
^200^ 
^00^ 

^Aiidorz^^^^ 1^200^ 
^elassa^^^Dhim ^20^) 
(^20^^ 

^o^ar^^aL (2^^) 
(^20^) 
(^2008) 

^ana^^^^ (^200^ 
(^200^ 
^200^ 

^^^FTi^^enko ^200^ 
(^200^ 

(^20^) 

^ucchiar^^^^^ (^20^^) 
Martin-Carrillo et al. (2014) 


Note. — a: For type I and II, ai is the decaying slope of the reverse shock emission. For others, it represents the rising slope of the forward 
shock emission, b: 0:2 is the decaying slope of the forward shock emission, c: Peak time in unit of s. For type I and II, it is the peak time 
of reverse shock emission. For others, it is the peak time of forward shock emission, d: Flux at peak time in unit of erg cm“^ s~^. 
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Fig. 3. — Optical lightcurves and fitting results for different types in selected sample, (a) Type I and II; (b)-(d) Type III; (e) Type IV; (f) 
Type III/IV. The majority of the data are collected in terms of observed magnitudes. Since most data are in the R band, we first calibrate 
the data from other wavelengths (“X” band) to the R band with the expression m/^ = mx — 2.5/3o log]^Q(Ai?/Ax) + 2.51ogiQ(/o ,R/f0,x), 
where /Sq is the the optical spectral indices (assuming Fb oc v~^o being satisfied in optical band), and /o is the absolute spectral irradiance 
for m = 0.0 within relevant magnitude system. An optical spectral index /3 q = 0.75 is adopted when 0q is not available JWang et al.ll201^ . 
1201511 . We then convert the R band magnitudes to the flux in units of erg cm~^ s~^ with the expression Fji = 

where Xr is the mean wavelength in R band. Galactic extinction correction is made to the data by using a reddening map presented by 
ISchlegel et al.l ll 199811 . 
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mation history. We ad opt a parame t erized GRB 
rate model proposed bv lYiiksel et al.l (|2008l l. 


Rgrb = Po 


(1 + ^) 


-34 


l + zV 
5000 ) 


+ 


1 + z 


35 


-I - 0.1 


The number of GRBs occurring per unit (observed) 
time in a comoving volume element dV{z)/dz is 
then 

dN Rgrb{z) dV{z) 


dtdz 


1 + z 


dz 


( 8 ) 


where the (1 + ^) factor accounts for the cosmolog¬ 
ical time dilation, and dV(z) /dz is given by 

dV(^) ^ _c _ IttDI _ 

dz Ho{i + zf[nM{i + zf + nA\y^' 
for a flat AGDM universe. 


• We assume that the electron spectral index p is 
the same for reverse and forward shock. Recent 
investigations suggest that the distribution of p is 
likely a Gaussian distribution, ranging from 2 to 
3.5, with a typica l value 2.5 fe.g. lLiang et al.l[2013l : 
IWang et al.ll2015ll . To test the influence of p value 
on the final results, the distribution of p is taken to 
be a Gaussian distribution with standard deviation 
0.2. Three values are tested for the mean value (p), 
i.e., 2.3, 2.5, 2.7. 

• The distribution function for the number den¬ 
sity of the ISM medium is still with large un- 
certainty, roughly ranging f r om 0.1 to 100 cm“^ 
(iPanaitescu fc Kumad l200ll 1200211 . Here we as¬ 
sume a Gaussian distribution in log space for the 
ISM density. The standard deviation is fixed 
as 10° ® (four orders of magnitude coverage for 
3a) and three mean values (h) are tested, i.e., 
1 , 10, lOOcm-3. 

• We assume that the fractions of shock energy that 
go to electrons (el’d) are the same for the reverse 
and forward shocks. The distribution of is 
poorly constrained in the literature, but due to the 
energetics consideration, in the past, a convention 
value 0.1 has been assumed in most studies. Here 
we assume a Gaussian distribution in log space 
for e^d with 10°'^ being the standard deviation^. 
Three values are tested for the mean value (^’■^), 
i.e., 0.001, 0.01, 0.1. 

• The fractions of shock energy that go into magnetic 
fields (e^'^) are very likely different, since the mag¬ 
netization degree of the ejecta material tends to be 
larger than in the ISM medium. We define'^ 


It has long been suggested that has a very 
wide distribution, ranging from 10~^ to 10~ ^ 
(|Panaitescu fc Kumail 12001 : ISantana et al.l [20T3l . 

In the simulations, Cg is generated through a uni¬ 
form distribution in log space with four ranges 
(each covering 3 orders of magnitude), i.e., 10“^ ^ 
10-1, 10-5 - 10-2, 10-5 - 10-5, lo-i" - 10-1. 

TT-s is generated through a Gaussian distribution 
with mean values (TZb) 1, 10 , 100, and standard 
deviation 1, y/TO, VlOO respectively. The value of 

f 

eg could be calculated with Cg and TZb straight¬ 
forwardly. 

• Gonsidering the power-law pr operty of the lum i- 
nosity function for GRBs le.g. iLiansr et al.l[2007|l . 
the kinetic energy of the GRB ejecta is generated 
with a power-law distribution. The minimum and 
maximum energy valu e are fixed as lO^i * erg and 
lO^'i erg respectively (|Zhang et al.ll200^ . Three 
power-law index ge are tested, i.e., 0.2,0.5,1. 

• The initial Lorentz factor of the GRB ejecta is also 
sugg ested to have a w ide range, from 50 to 500 
le.g. iLiang et al.ir2013ll . Here we generate Tq with 
a uniform distribution in the log space. We test 
three combinations of the minimum and maximum 
values, i.e., 50 ^ 300, 100 ^ 500 and 50 ^ 500. 

• The observed shell width essentially shares the 
same distribution with the GRB duration, which 
could be well described wit h a Gaussian di stribu- 
tion in the log space^ (e.g. iQin et alllMll . The 
observed shell width is generated with a Gaussian 
distribution in log space with standard deviation 
10°'5. We test two mean values (Aq), i.e., lO^^ cm 
and 10^2 cm. 

3.2. Simulation Results 

Given a set of distribution functions for each afterglow 
parameter, we run Monte Carlo simulation for 10000 
times 5, and we analyze the thus obtained distributions 
of fractional ratios between different types of lightcurves. 
In Figures |31[S] we plot the simulation results for selected 
situations which are relevant for illustrating the main 
conclusions, which can be summarized as follows: 

• The fraction ratios between different lightcurve 
types depend sensitively on two key parameters, 
e^d and TZb- The value of 7 ?._b characterizes the 
balance between reverse shock dominated cases (I 
and H) and forward shock dominated cases (HI and 
IV). Increasing TZb can significantly increase the 
proportion of Type I and H. The value of e^d es¬ 
sentially determines the internal coordination be¬ 
tween Type I and H, or Type HI and IV. Smaller 
e^d gives more Type HI and Type H. 


TZB^e'Blei- ( 10 ) 

^ We also tested larger standard deviation values. It turns out 

r f 

that if the standard deviation for is too large, the observational 
results, especially the internal coordination between Type I and II 
could never be reproduced. 

^ This definition is different from the original definition of Zhang 
et al. (2003), who defined 'TZb = BrjBf, which is the square root 
of the 'TZb defined in this paper. 


• When W^d — q.I, as shown in Figure |4^-l4]:, TZb 
should be larger than 10 but smaller than 100, oth¬ 
erwise the proportion of Type I and H is either too 

® We only consider long GRBs here, since the collected sample 
are essentially all long GRBs. 

® The number of runs for each simulation is determined by bal¬ 
ancing the computation time consumption and the resulting con¬ 
vergence. 
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£g distribution range Eg distribution range Eg distribution range 

(g) 4 = 0.1,TIb = 10, p= 2.7 (h) = 0 . 1,113 = 10, Tq = 100 ~ 500 (i) e{ = 0.1,Tes = 10 ,Tq = 50 ~ 500 


Fig. 4.— Stacked fractions of different lightcurve types for observational data and selected simulation results. Each panel corresponds 
to one specific simulation setup. Unless specified under the subfigure, the general setup values for different parameters are as follows (see 
details in section 3.1): n = l,p = 2.3, Fq = 50 ~ 300,0b = —0.5, Aq = 10*^^. Five stacked histograms are presented in each panel. The 
first four histograms are for different ranges (left to right: 10“^ ~ 10“'^, 10“® ~ 10“^, 10~® ~ 10“^and 10“^ ~ 10““^) and the last 
one represents the observational results. 


small or too large to reproduce the observational 
data. On the other hand, the observed fraction of 
Type II is much larger than Type I, which is in 
contrast with the simulation results. Varying the 
value of TLb does not help to adjust the fraction 
ratio between Type I and II. 

• Keeping e^’-^ of order of 0.1 and TZb oi order of 10, 
we also checked if the observational results could be 
reproduced by varying other parameters. Since the 
simulation results depend sensitively on the value 
of and TZb, for better testing the effects of 
other parameters, we fix the value of as 0.1, the 
value of TZb as 10, when the distribution functions 
of other parameters are being varied. The mean 
value of number density n is varied from 1 to 10 
and 100; the mean value of electron index p is var¬ 
ied from 2.3 to 2.5 and 2.7; the distribution range 


of initial Lorentz factor Tg is varied from 50 ^ 300 
to 100 ^ 500 and 50 ^ 500; the power law index of 
kinetic energy distribution function is varied from 
0.5 to 0.2 and 1; and the mean value of the initial 
shell width Ag is varied from 10^^ cm to 10^^ cm. 
As shown in Figure varying the distributions 

of these parameters does not affect the results too 
much and hence does not help to solve the incon¬ 
sistency of the ratio between Type I and Type II. 

• Fixing the distributions for all other parameter, the 
observations can be easily reproduced as long as 
the value of 4’^ reduced by one order of mag¬ 
nitude i.e., 4’^ = 0-01 (when IZb = 100 as shown 
in Figure [Sf). The constraint on Cg is not strong, 
but smaller values of ranging from 10“® to 10“^ 
seems to be more favorable. 
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£g distribution range Eg distribution range Eg distribution range 

(a) = 0.1, Tea = 10, as = -0.2 (b) el = 0 . 1,113 = 10, ob = -1 (c) el = 0 . 1,113 = 10, Aq = lO^^ 



Eg distribution range Eg distribution range Eg distribution range 

(d) = 0.01,:^s = 1 (e) = 0.01,7 ?.b = 10 (f) = 0.01,= 100 



Eg distribution range Eg distribution range Eg distribution range 

(g) = 0.001,7 ?.b = 1 (h) = 0.001,-^B = 10 (i) = 0.001,-^B = 100 


Fig. 5.— Same a Figure!^ but for different parameters. 


• When e^’-^ = 0.001, although the fraction of Type I 
and II could be consistent with the observations as 
long BuS 1 Zb is large enough, the fraction of Type IV 
is too small (or even completely disappear), which 
is inconsistent with the observations. 

In summary, the simulation results indicate that our 
morphological analysis for early optical afterglow is able 
to efficiently constrain the microscopic parameters, e.g., 
eg and TZb- To reproduce the current observations, 
= 0.01, IZb = 100 and relatively smaller values of 
Cg is favored, which can be understood as follows: in 
the observational data, the fraction of Type II is larger 
than Type I, inferring that the peak of the forward shock 
emission is easily suppressed by the reverse shock com¬ 
ponent. On the other hand, the fraction of Type III is 
larger than Type IV, even when all the bursts of over¬ 
lap type belong to Type IV. As illustrated in Figure [H 
both these items of observational evidence can be ex¬ 
plained if the forward shock component is in the FS II 
case {twilitx) < t'opt), which favors a smaller value of e{. 


If el becomes smaller, the forward shock emission in the 
optical band becomes stronger, so that a larger value of 
TZb and a relatively smaller value of eg is required to 
maintain the balance between reverse shock dominated 
cases and forward shock dominated cases. 

4. DISCUSSION 

A practical scheme of morphological analysis for GRB 
early optical afterglows and its ability to constrain af¬ 
terglow parameters has been illustrated in the last two 
sections. We have applied this method to the currently 
available observational results, and have derived con¬ 
straints on the relevant microscopic parameters. In the 
following, we will discuss some of the challenges facing 
this method and the caveats on our constraint results. 

The greatest challenge for the morphological analysis 
method arises from the sample selection. It is difficult 
to achieve the completeness of a certain sample, unless 
a sufficiently large number of triggered GRBs can be 
rapidly followed-up in the optical band. On the other 
hand, systematic uncertainties could become large, and 
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would be difficult to remove, if the afterglow follow-ups 
are obtained through different telescopes. A future ded¬ 
icated facility with rapid response ability and wide field 
of view could help with these issues, and this is a key ele¬ 
ment in the Chinese-French missi on SVOM, the G round 
Wide Angle Cameras (GWACs) (jPaul et al'll201111 . 

Another challenge comes from the process of assign¬ 
ing the observed lightcurves into relevant categories. To 
better identify the Type III and Type IV, sufficient data 
points in the rising phase are required, while to precisely 
distinguish Type I from Type II, observations in the de¬ 
caying phase need to be dense enough. Multi-color ob¬ 
servations during the follow-up phase are essential to ad¬ 
dress this challenge. 

The theoretical scheme for determining lightcurve cat¬ 
egories is based on the standard synchrotron external 
shock model. Despite its great success, the standard 
model has some limitations that sometimes hinder a pre¬ 
cise description of GRB afterglows. For instance, the real 
evolution of may deviates from a power-law behavior 
when t is around tx, so that both the reverse shock and 
forward shock lightcurves should have a smooth transi¬ 
tion around the peak, especially when equal arrival time 
effects are considered. These deviations may affect our 
results over some limited range of parameter spaces, e.g., 
when Pf (</) is close to {t^). Such effects may aver¬ 
age themselves out, as long as the simulated sample is 
large enough. In principle, one can use numerical sim¬ 
ulations to calculate more precise lightcurves for given 
set of parameters, but this will dramatically increase the 
computation time while most of the calculations are re¬ 
dundant for the purpose of morphological analysis. 

Due to the limitations of the current facilities, the sam¬ 
ple selected in this work is still incomplete in some sense. 
As mentioned in section 2.3, only II4 swift bursts have 
optical follow-up within 500 s, and some of them are 
hard to classify because of lack of sufficient data points 
in the rising phase. The incompleteness may cause some 
uncertainty in the parameter constraint results, but the 
general tendency of our results should be reliable in or¬ 
der of magnitude, e.g., should be in order of O.OI and 
TZb should be in order of 100. 

The analysis in this work is designed for a homogeneous 
interstellar mediu m. For GRBs occurrin g in a wind type 
environment le.g. iGhevalier fc LilllQQOl l. the lightcurves 
are easy to distinguish from what is discussed here, and 
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